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Strong lensing time delay cosmography has excellent complementarity with other dark energy 
probes, and will soon have abundant systems detected. We investigate two issues in the imaging 
and spectroscopic followup required to obtain the time delay distance. The first is optimization of 
spectroscopic resources. We develop a code to optimize the cosmological leverage under the con¬ 
straint of constant spectroscopic time, and find that sculpting the lens system redshift distribution 
can deliver a 40% improvement in dark energy figure of merit. The second is the role of systematics, 
correlated between different quantities of a given system or model errors common to all systems. 
We show how the levels of different systematics affect the cosmological parameter estimation, and 
derive guidance for the fraction of double image vs quad image systems to follow as a function of 
differing systematics between them. 


I. INTRODUCTION 

Cosmographic, or geometric, methods such as 
distance-redshift relations provide key insights into the 
nature of our universe. The Type la supernova luminos¬ 
ity distance-redshift relation revealed that the cosmic ex¬ 
pansion is accelerating [1, 2], with the physical cause de¬ 
noted as dark energy. The cosmic microwave background 
radiation anisotropies and baryon acoustic oscillations in 
galaxy clustering are other probes that have at least a 
substantial geometric component (modulo dark energy 
perturbations or coupling to matter). Cosmic redshift 
drift has been recognized since the 1960s as a potential 
cosmographic probe, though not yet measured [3-6] . The 
strong gravitational lensing time delay distance-redshift 
relation was also proposed in the 1960s [7] and in the 
last few years has matured to the stage of being used as 
a cosmological probe [8, 9]. 

The development of strong lensing distances is a par¬ 
ticularly interesting advance since unlike the standard 
distance-redshift relations the measured time delay is a 
dimensionful quantity, and the time delay distance is 
comprised of the ratio of three distances. This makes 
it sensitive to the Hubble constant Hq. Also because of 
this ratio the time delay distance has an unusual depen¬ 
dence on dark energy properties and has high comple¬ 
mentarity with the usual distance probes [10, 11]. On¬ 
going and future surveys such as from the Dark Energy 
Survey, Large Synoptic Survey Telescope, and Euclid and 
WFIRST satellites have incorporated strong lensing time 
delays into their suite of cosmological probes. 

Here we examine two aspects of implementation of time 
delay distances into such surveys, focused on trades and 
optimization of the followup resources required to obtain 
a robust distance-redshift relation. In particular the wide 
field imaging surveys must be supplemented with spec¬ 
troscopy to obtain accurate redshifts of lens and source, 
and to constrain the lens mass model. Since spectroscopy 
is time intensive, and not part of some of the planned 
surveys, we consider how to efficiently allocate the addi¬ 
tional resources among the large numbers (1000-10000) 


of strong lens systems that will be found. 

In Sec. H we review the basics of strong lensing time 
delays and the types of observations necessary to mea¬ 
sure the time delay distance-redshift relation. We de¬ 
velop in Sec. HI an optimization procedure for the cos¬ 
mological leverage of the data under the constraint of 
fixed resources such as total spectroscopic time. This 
addresses questions of followup of low vs high redshift 
systems. The influence of systematic uncertainties is in¬ 
vestigated in Sec. IV, along with questions such as how 
to trade between different populations of lens systems, 
such as ones with double images vs quad images. We 
summarize and conclude in Sec. V. 


II. MEASURING TIME DELAY DISTANCES 

The time delay distance can be thought of as the focal 
length of the lensing, and depends on the distances be¬ 
tween observer and lens Di , observer and source Dg , and 
lens and source Du . The time delay between two images 
of the source comes from the geometric path difference of 
the light propagation and from the differing gravitational 
potentials experienced. In summary, the time delay dis¬ 
tance is 
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where zi is the lens redshift. At is the observed time delay, 
and A(j> is the Fermat potential difference modeled from 
the observations such as image positions, fluxes, surface 
brightness, etc. 

For strong lensing time delay cosmography, the source 
should be a bright, time varying object such as an active 
galactic nucleus (AGN) and the lens is generally a fore¬ 
ground galaxy (as cluster lenses are harder to model). 
See [12-22] for further details on strong lensing time de¬ 
lays as a cosmological probe. 

Wide field surveys such as from the Dark Energy Sur¬ 
vey, Large Synoptic Survey Telescope, and Euclid and 
WFIRST satellites will be superb tools for finding large 
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samples of strong lens systems. Due to their repeat ob¬ 
servations, they can also monitor the image fluxes over 
several years to measure the time delay At. This may be 
supplemented with further cadenced observations from 
external programs, along the lines of the highly success¬ 
ful COSMOGRAIL program [23]. 

The Fermat potential is constrained by the rich 
data of the images, but this works best with additional 
high resolution imaging, currently supplied by the Hub¬ 
ble Space Telescope, and for future surveys possibly by 
the James Webb Space Telescope (JWST) and ground 
based adaptive optics. For the lens part of the Fermat 
potential, the lens mass modeling requires constraint by 
measurement of the galaxy velocity dispersion through 
spectroscopy. This also plays a key role in breaking the 
mass sheet degeneracy [9, 24, 25]. Similarly, spectroscopy 
obtains the redshifts of lens and source. (The velocity 
dispersion is also crucial for the possibility of using time 
delay lensing to obtain the usual angular diameter dis¬ 
tance [21, 26].) 

Thus these essential followup resources must be sought 
in order to derive the strong lensing cosmological con¬ 
straints from the wide field imaging survey. Since these 
are generally external to the wide field survey, and require 
application for highly subscribed telescope time, they can 
become a limiting factor in the science return. We con¬ 
sider here the optimization of cosmological leverage given 
a finite followup resource. In the next section we present 
calculations specifically dealing with spectroscopy, but 
the optimization procedure is quite general. 

Since any one cosmological probe has particular de¬ 
generacies between parameters, we combine the strong 
lensing distances with cosmic microwave background and 
supernova distances. Strong lensing was shown to have 
great complementarity with these probes [11], and these 
data will exist at the time of the wide field surveys (in¬ 
deed supernovae as distance probes are another compo¬ 
nent of the surveys). We adopt a Planck quality con¬ 
straint on the distance to last scattering (0.2%) and phys¬ 
ical matter density (0.9%). For supernovae we use 

a sample of the quality expected from ground based sur¬ 
veys: 150 supernovae at z < 0.1, 900 from z = 0.1-1, 
and 42 from z = 1-1.7, with a statistical uncertainty of 
0.15 mag and a systematic of 0.02(1 -|- z) mag added in 
quadrature to each 0.1 width bin in redshift. 

We perform a Fisher information analysis to estimate 
the cosmological parameters of the matter density Dm, 
dark energy equation of state present value wq and a mea¬ 
sure of its time variation Wa, reduced Hubble constant h, 
and a nuisance parameter A4 for the supernova absolute 
magnitude. The hducial cosmology is flat ACDM with 
Dm = 0.3, h = 0.7. 


III. OPTIMIZING SPECTROSCOPIC 
FOLLOWUP 

Since spectroscopic time is restricted, and generally 
requires arrangements outside the main survey, it is ad¬ 
vantageous to treat it as a limited resource and optimize 
its use. We consider it as a fixed quantity, and seek to 
maximize the cosmological leverage of the measured time 
delay distance given this constraint. To do so, we exam¬ 
ine the impact of sculpting the redshift distribution of the 
lenses to be followed up. The spectroscopic time is domi¬ 
nated by measurement of the lens galaxy velocity disper¬ 
sion. However our methodology is general and similar 
results should occur for any measurement to some given 
signal to noise. For example, one might instead optimize 
a resource capable of high resolution imaging, to map the 
distorted source images, such as JWST or ground based 
adaptive optics time. The principles are the same. 

For specificity, we concentrate on fixed spectroscopic 
time for the sample of lenses. To measure a redshift, 
or the galaxy velocity dispersion, requires good signal to 
noise data of line fluxes. Consider the following illustra¬ 
tive calculation. The signal scales with the number of 
photons from the spectral feature to be measured, hence 
the fluence times the exposure time. We will take the 
contribution of other sources of photons to be dominated 
by a redshift independent contribution, times the expo¬ 
sure time, so the noise, i.e. the fluctuations, goes as the 
square root of exposure time. This gives 

N ^ 

We emphasize that we are presenting an illustrative 
methodology: true survey optimization will depend on 
many sources of noise and the survey specifications such 
as instrumental properties, scanning strategy, etc. This 
requires a full exposure time calculator and is beyond 
the scope of this article, but we will see below that our 
simple, heuristic approach matches some known results. 

Since the fluence is just the flux divided by the photon 
energy, we lose one less factor of 1 -I- z than the usual 
flux-redshift relation, i.e. T oc {1 -\- z )/where is 
the luminosity distance. So to achieve a desired constant 
signal to noise threshold requires the exposure time to 
vary as 

1 

texp (Tj^ 

where = {1 + z)~‘^dL is the angular diameter distance. 
At redshifts z ~ 1-2, the angular diameter distance in 
our universe stays close to constant, and we recover the 
result [27] that spectroscopic exposure time becomes in¬ 
creasingly expensive with redshift as roughly (1 -|- z)®. 

At lower redshift, the slope is steeper. We will be in¬ 
terested in a range around z ~ 0.5. However, as expo¬ 
sure time gets smaller, other noise contributions enter as 
well as overheads such as telescope slewing and detector 
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readout time. Therefore we adopt a reasonable approx¬ 
imation that the spectroscopy cost goes as t oc (1 -I- 2 )’' 
with r = 8. We have checked that using instead Eq. (3) 
plus a constant overhead makes no significant difference 
in our results. We emphasize again that the key point is 
that for any texp(-z) from an instrument’s exposure time 
calculator, the optimization code described will produce 
results under a fixed resource constraint. 

The next step, given the resource constraint, is to 
choose the quantity to optimize. We take this to be 
the dark energy figure of merit (FOM), the area of a 
confidence contour in the dark energy equation of state 
plane, marginalized over all other parameters. The dark 
energy equation of state w(a) = wq + 'M^a(l — o) fits 
a broad range of models, and is accurate in recreat¬ 
ing distances to the 0.1% level [28, 29]. The parameter 
wq measures the present equation of state and Wa the 
time variation, with a = 1/{1 + z) the scale factor, and 
FOM = (det COV[wo, In the next section we 

also consider the effect of optimization on the determi¬ 
nation of the Hubble constant. 

At different lens redshifts, the time delay distance 
has different sensitivities to the cosmological parameters. 
(This is true for the source redshift as well, but we fix 
Zs = Szi for simplicity; reasonable variations of this ratio 
have little effect on the cosmological sensitivity [11].) So 
the question is whether the extra expense of spectroscopy 
of higher redshift lenses overcomes their possibly greater 
leverage. 

To optimize the redshift distribution we begin with a 
uniform distribution in lens redshift (recall we are most 
interested in spectroscopy of the lens galaxy to obtain 
its velocity dispersion, used to constrain the lens mass 
model). We take 25 time delay systems of 5% precision 
in each bin of redshift width dz = 0.1 over the range 
z = 0.1“0.7, for a total of 150 systems. This carries 
with it a certain total spectroscopic time, and that is the 
fixed resource constraint under which the optimization 
proceeds. 

This initial uniform distribution is perturbed by one 
system in each bin, one at a time, and resulting FOM is 
calculated. Each redshift also has a different time bur¬ 
den, and we compute the quantity 

, FOM(perturbed) - FOM , , 

for each redshift bin i, where Ati is the spectroscopic time 
required for a system at that redshift. The bin with the 
lowest FOMi, i.e. the least change in cosmological lever¬ 
age, has one system subtracted from it. The time saved 
is then reallocated to the other redshift bins, increasing 
the number of systems in every other bin j , weighted by 
FOMj. That is. 
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/iff -2^sub 

V 1 ff 


FOMj 


( 5 ) 


where sub is the bin from which Ansub = 1 systems are 


removed. This formula conserves the resource, i.e. spec¬ 
troscopic time. 

The FOM for the new distribution is computed, and 
the process iterates. The new distribution is perturbed, 
and again one system from the lowest leverage bin is re¬ 
moved (if this would cause the number in that bin to 
go negative, we use the next lowest leverage bin) and 
its time burden is reallocated. The iteration continues 
until convergence. As a final step we round the num¬ 
bers in each bin to the nearest integer, but this has less 
than 0.3% impact on the cosmology parameter estima¬ 
tion. This optimization method is computationally fast 
and efficient, and widely applicable to many astrophysi- 
cal studies with constrained resources. 

Figure 1 illustrates the results. The optimization in¬ 
creases the FOM by almost 40%, while keeping the spec¬ 
troscopic time fixed. The optimized redshift distribu¬ 
tion has a number of interesting properties: it is heav¬ 
ily weighted toward low redshift, with a single peak at 
higher, but not maximal redshift. Low redshift gives a 
decreased time burden, and still good cosmological lever¬ 
age, especially on the Hubble constant h, but also the 
dark energy parameters since the source redshift extends 
the cosmological lever arm [11]. To break covariances 
between parameters, the higher redshift bin is needed, 
but note it does not seek to maximize the range by tak¬ 
ing the highest bin since this has the greatest time bur¬ 
den. (Note that the pioneering cosmological optimization 
of [30] fixed the number of supernovae, not the observ¬ 
ing time, to find a peak at the redshift maximum. In 
[31] both constant number resource and constant spec¬ 
troscopic time resource with n = 6 were studied for the 
supernova distance probe.) The intermediate redshift 
bins, and the two highest redshift bins, at z = 0.5-0.6 
and 0.6-0.7, are zeroed out in the optimization. 

The FOM becomes a quite respectable 192, with deter¬ 
mination of ftm to 0.0035, Wo to 0.061, Wa to 0.22, and 
h to 0.0030. Each parameter estimation, as well as the 
FOM, is better than for either the uniform redshift distri¬ 
bution or a magnitude limited distribution derived from 
[14] using cuts on image and lens flux and image sep¬ 
aration (P. Marshall, S. Suyu private communication). 
We have tested reducing the redshift range to 0.1-0.4 or 
0.1-0.5, and obtain the same optimized distribution, i.e. 
the optimum really has only the low and mid redshift 
bin. We also shifted the range to z = 0.2-0.7 and found 
the optimized FOM dropped significantly, to 135, demon¬ 
strating the lowest bin is critical for dark energy as well 
as the Hubble constant. 


IV. INFLUENCE OF SYSTEMATICS 

Any cosmological probe must deal with systematic un¬ 
certainties, especially for next generation surveys where 
abundant numbers of objects drive down the statistical 
uncertainty. We investigate two of the manifold aspects 
of the impact of systematics. Again, a detailed treatment 
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FIG. 1. Histograms of the lens redshift distribution are shown 
for three cases: uniform, magnitnde limited, and optimized 
cosmological leverage distributions. All three have the same 
fixed resource constraint on the total spectroscopic followup 
time; each is labeled with the resulting dark energy figure of 
merit. 

would need to delve deep into survey and instrumenta¬ 
tion properties, and is beyond the present scope. 

A. Redshift distribution revisited 

Let us explore the effect of systematic uncertainties 
on the optimization carried out in Sec. III. The num¬ 
ber of lens systems in the lowest redshift bin approached 
300, which for individual system precision of 5% implies 
a required control of systematic bias at the 0.3% level. 
We do not yet know enough from current strong lens- 
ing observations and studies to know what is a realistic 
level, but strong efforts and advances in understanding 
are underway. For example, the blind data Time De¬ 
lay Challenge (http://timedelaychallenge.org) has 
already achieved 0.1-0.2% control of time delay estima¬ 
tion [32-35]. 

Therefore we study the impact of various levels of sys¬ 
tematic on the optimized redshift distribution and the 
resulting cosmological parameter estimation. We imple¬ 
ment the systematic as a floor, added in quadrature to 
the statistical uncertainty. 


where rii is the number in redshift bin i and each bin 
is treated independently. This model is commonly used 


in supernova distance-redshift relation studies [36] . (But 
see the next subsection for an alternative approach to 
systematics.) 

Figure 2 shows the optimized distributions, subject to 
the resource constraint, for different levels of systematic. 
As the systematic level increases, it is less advantageous 
to put a large number of systems in a given redshift 
bin and additional systems diffuse into neighboring bins. 
This fills in the intermediate redshift gap and also pushes 
some lens systems to high redshift, making the distribu¬ 
tion closer to uniform. 



z 


FIG. 2. The resource-constrained optimization of the lens 
redshift distribution is shown as a function of coherent dis¬ 
tance systematic floor. This floor prevents large numbers of 
systems in a redshift bin from improving the accuracy. 

Both the presence of the systematic and the redistri¬ 
bution of the data away from the zero-systematic opti¬ 
mum lowers the FOM. Figure 3 plots the FOM vs the 
systematic level. At 0.2% systematic, the FOM has de¬ 
creased by only 8% relative to zero systematic case, but 
larger systematics impact the cosmology more severely. 
For high enough levels, there is little difference in lever¬ 
age between the optimized and uniform (or magnitude 
limited) distributions. 

To make sure that optimizing for dark energy FOM 
also helps improve other cosmological parameters, we 
show in Fig. 4 the constraint on the Hubble constant h. 
For zero systematic the optimization actually improves 
the estimation by a factor 2 (not just 40% as for the 
FOM). This is because the increased low redshift sample 
is particularly useful for the Hubble constant. Even for 
higher levels of systematics the optimization continues to 
give added leverage on h. 
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FIG. 3. Dark energy figure of merit is plotted for the opti¬ 
mized cosmological leverage and uniform lens redshift distri¬ 
butions as a function of systematic floor. For high systemat- 
ics, the optimized distribution has little extra leverage, but 
at low systematics the improvement can approach 40%. 



FIG. 4. As Fig. 3 but for the Hubble constant constraint 
(note that here low values are better). At high systematics, 
the optimization still improves the results by almost 10%, 
while at low systematics the gain is a factor 2. 


B. Model systematics 

Statistical uncertainties in the time delay distance arise 
from measurement imprecision, but the measurement er¬ 
rors can have systematic components as well. These take 
two main forms: correlated errors between eleirrents in a 
given lensed system, and irrodel errors that are common 
between different lens systems. These will give diago¬ 
nal and offdiagonal contributions to the distance error 
matrix. Below we give an illustrative exhibition of the 
effect of such systematics; again, actual survey analysis 
requires a more detailed and sophisticated treatment but 
this demonstrates the main points. 

Let us write the time delay distance as 

^’At = -DAt(At, A^(m, 0, At, V, z), Kext) , (7) 

where the Fermat potential depends on image magni¬ 
tudes m, positions 6, time delays At, lens and image 
redshifts z, and lens velocity dispersion v. Additional 
mass along the line of sight affects the modeling through 
the external convergence Kext- 

Ref. [9] demonstrates that over the angular range im¬ 
portant for the images, i.e. near the Einstein radius of 
the lens, the Fermat potential scales with the projected 
lens mass profile slope 7 '. Even when the galaxy lens 
profile has a multicomponent composition [37] of a Hern- 
quist stellar core plus a Navarro-Frenk-White dark mat¬ 
ter profile, over the angular range of interest the slope 7 ' 
captures the profile dependence. Uncertainties or mises- 
timation in 7 ' then lead to dispersion or systematics in 
the Fermat potential and hence the distances. There is 
also a mass sheet degeneracy due to mass along the line 
of sight. Both these effects can be incorporated through 

[9] 

())Ri <^(1 - Kext)(7' - 1) , (8) 

in the vicinity of the standard profile slope 7 ' = 2. Note 
that systems measured to date have rms dispersion of 
< 5% in 7 ' so this form is accurate. 

We can now write the time delay distance as 


(j){m, 9, At, V, z) (I - «;ext)(7' - 1) 

The error propagation to the time delay distance, for 
accurate image flux and position measurements, is then 

5D = DtSt + Dy6v + Dz5z + Dk.Sk + , ( 10 ) 

where we use the simplifying notation D = D/^t, 
Dx=dD/dx, and 7 = 7 ', k = Kext, t = At. 

A reasonable first approach, based on current data, is 
that the error budget will be dominated by the time delay 
estimation, external convergence, and lens mass profile 
and velocity dispersion. In this case the diagonal entries 
in the error matrix (i.e. for a single lens system) would 
be 

Cdd ^ + Dial + Dial+ Dy^ 

+DyD^{6v Sj) + DkDj{SkSj) . ( 11 ) 
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This reflects the individual errors, plus the correlated 
errors between v and 7 in the lens density profile (e.g. 
whether v is measured at the appropriate place in the 
profile), and between 7 and k in the profile-mass sheet 
degeneracy. We emphasize this is illustrative. 

The model errors enter in the off diagonal elements. 
Recall these correlate two different lens systems at dis¬ 
tances D and D'. Here the error matrix gets the contri¬ 
bution 


Cdd' = DtD[{5t 5t') + DyD'^{dv Sv') + DkD'^{Sk6k') 

+ DjD'^{S'yS'j') + DyD'^{Sv5^) + D'^D^{5v S'y) 
+ D^D'^{SkSj') + D'^D^{Sk' S-f) . ( 12 ) 

Several of the derivatives can be written in a straightfor¬ 
ward manner using Eq. (9): 


^ D ^ D 

Dt — — ; — — 


At 


D,, = 
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1 ^ext 


7-1 

dd> 


(13) 
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The approximation sign in Dy represents the result for a 
singular isothermal sphere. 

Once we are in the regime of thousands of strong lenses 
[14, 15], the \fN statistical reduction will be dominated 
by the residual systematics. It is crucial to identify the 
effects of these systematics on cosmological results, and 
where the greatest leverage lies in controlling them. This 
has significant interplay with our previous analysis of 
constrained followup resources. That is, we want to iden¬ 
tify where best to concentrate the limited resources, e.g. 
on long time delay, quad image systems. One example is 
that double images often occur at different radii where 
the lens slope profile may vary, while quads suffer less 
from such a systematic. Quad systems may have better 
precision from the extra measurement constraints, but 
possibly also an increased opportunity for differential mi- 
crolensing or varying external convergence. 

For a tractable first approach, we note that Eqs. (11), 
(12), and (13) have the property that the uncertainties 
often enter into the error covariance matrix as logarith¬ 
mic fractional quantities, i.e. 5^'/— 1) or 5k,/{1 — k), 
times the time delay distance. We adopt the Ansatz that 
these fractional systematic uncertainties have a scaling 
with redshift s = s* [(1 -I- z)/(l -I- 2 *)]", where z* = 0.4 is 
the midpoint of the redshift range, and a possible pop¬ 
ulation dependence, where the errors in double systems 
may differ from those in quad systems. Since the frac¬ 
tion of doubles vs quads detected in a survey changes 
with redshift, this causes a population drift in a manner 
similar to supernova subtype evolution. 

The overall covariance matrix will then have two en¬ 
tries in each redshift bin, for doubles and quads, with 
distinct statistical and systematic errors. We sum up all 
the statistical errors from the time delay estimation, mass 
profile slope, etc. to give diagonal entries of 
in the ith lens redshift bin. Here Di is the time delay dis¬ 
tance of Eq. (1). The systematic errors S[d,q}{zi) Di also 


contribute to the diagonal elements and moreover their 
correlations produce offdiagonal entries. While in actual 
data analysis one might not bin the data, and the er¬ 
ror model will become more sophisticated over time, this 
approach using five parameters (crd, aq, Sd, Sq, and n if 
desired) is tractable and gives important first indications 
of the effect of systematics. 

We first discuss direct redshift evolution of the system¬ 
atics, and then the influence of population drift. 

Figure 5 shows the influence of the redshift dependence 
of the systematic. In order to focus on the redshift evo¬ 
lution of the systematic, we here take the uniform red¬ 
shift distribution of Sec. HI (despite its lower FOM), and 
treat all populations (doubles or quads) as having 5% 
statistical precision in distance and a fractional distance 
systematic s per redshift bin. The evolution in s with 
redshift could arise from, e.g., decreased signal to noise, 
and hence more uncertain modeling, of higher redshift 
lens systems (and a longer path length so greater uncer¬ 
tainty in the projected mass along the line of sight). 



Systematic s. (%) 


FIG. 5. The dark energy figure of merit (solid curves) and 
Hubble constant uncertainty (dashed curves) are shown as a 
function of systematic error amplitude s* at z* = 0.4, for 
three different redshift evolutions. The systematic scales as 
(1 -I- z)" for n = 0, 1, 2. Note that high figure of merit is 
good, and low (t(/i) is good. 


We see that the specific systematic redshift evolution 
model, as opposed to the mere presence of the systematic, 
has a modest effect, with the variation among n = 0 — 2 
affecting the dark energy figure of merit by ^ 20 % and 
determination of the Hubble constant by ^ 50% at the 
highest systematic levels. Note that the evolution does 
change the covariance between parameters; this is respon- 
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sible for the n = 1 and n = 2 curves crossing at high 
systematic. Although wq and Wa are both better deter¬ 
mined in the n = 1 case, the FOM is not, due to their 
altered covariance. As expected, a constant fractional 
systematic has more of an effect on /i, and so less on the 
FOM. 

As a next step, we include separate systematics for the 
populations of double and quad image systems. Again 
the motivation is that these have different levels of con¬ 
straints on the lens model from the observations, e.g. 
image positions, flux ratios, number of time delays. To 
focus on this population aspect, we keep the explicit sys¬ 
tematic amplitude independent of redshift (n = 0; recall 
we just saw that the cosmology constraints were fairly in¬ 
sensitive to n anyway) but incorporate population drift 
with redshift. 

The distribution of doubles and quads we use is the 
magnitude limited sample of Sec. Ill, based on [14] (again 
for illustration, despite its lower FOM). This arises from 
the constraints that 1) image separation is larger than 1", 
2 ) the quasar images are brighter than i-magnitude 20.8, 
and 3) the lens magnitude is brighter than 22. These 
help ensure that images and arcs can be well resolved, 
time delays can be accurately measured, and the lens ve¬ 
locity dispersion can be measured with reasonable use 
of resources. The ratio of doubles to quads varies from 
roughly 9 to 4 from the low to high end of the redshift 
range, with statistical scatter in numbers included. The 
drift in the proportion leads to an effective redshift evo¬ 
lution in the systematic uncertainty, analogously to how 
population drift of supernovae subtypes with slightly dif¬ 
ferent absolute magnitudes engenders supernova magni¬ 
tude evolution (see [38, 39] for detailed treatment of the 
propagation of this effect into cosmology constraints). 

Figure 6 quantifies how population drift between the 
double image and quad image systems propagates into 
cosmological constraints as the amplitude and ratio of the 
double and quad systematic errors varies. The top panel 
fixes Sq = 1% and varies s^, while the bottom panel fixes 
Sd = 4% and varies Sq. As the systematic level increases, 
the dark energy figure of merit decreases and the uncer¬ 
tainty in the determination of the Hubble constant in¬ 
creases, as expected. However, when the systematic error 
of the doubles exceeds their statistical uncertainty, then 
the information from the doubles saturates, the quads 
dominate the leverage, and the FOM and <j{h) level off 
(see top panel). This implies that systematics in doubles 
have a natural “knee” - defining when the limited ob¬ 
servational resources have greater leverage when used to 
accurately characterize the rarer quad systems. The bot¬ 
tom panel shows that as long as the quad systematic is 
lower than its statistical uncertainty, then more accurate 
measurement of quads leads to more stringent cosmolog¬ 
ical constraints. Such an analysis can inform the optimal 
use of resources for strong lensing time delay distances 
as a cosmological probe. 

Figure 7 depicts the results for simultaneous variation 
of the double and quad systematics, showing contours 



FIG. 6. [Top panel] The dark energy figure of merit (solid 
curve) and Hubble constant uncertainty (dashed curve) are 
shown as a function of fractional systematic error Sd in double 
image systems, for a fixed quad image systematic Sq = 0.01. 
[Bottom panel] As the top panel, but as a function of Sq for 
fixed Sd = 0.04. 


of constant FOM. At high systematics level, having the 
two contributions be comparable gives the best FOM. 
However at low systematics level we are below the quad 
statistical uncertainty. This gives two approaches for im¬ 
provement: either improving the doubles systematic or 
allocating more resources to followup more quad systems 







and bring down their statistical error. From the distance 
between the contours, we see that if low quad systematics 
can be achieved, then considerable improvement in dou¬ 
bles systematics is required for significant improvement; 
thus following up more quad systems seems a better op¬ 
tion in this case. 



s, (%) 


FIG. 7. Contours of constant FOM are plotted as a function 
of the systematic levels in double and quad leasing systems. 
The zero systematics point has FOM=126. 


V. CONCLUSIONS 

The strong gravitational leasing time delay distance- 
redshift relation is a geometric probe of cosmology. It 
has two particularly valuable characteristics: being di¬ 
mensionful and hence sensitive to the Hubble constant 
Hq, and being a triple distance ratio and hence with dif¬ 
ferent parameter degeneracies that make it highly com¬ 
plementary to other distance probes such as the cosmic 
microwave background or supernovae. Moreover, the ob¬ 
servations and modeling are rapidly advancing, enabling 
it to place cosmological constraints of significant lever¬ 
age, comparable to other methods. 

We considered the question of the followup resources 


needed to complement the forthcoming strong lensing 
imaging surveys that detect and monitor the lens sys¬ 
tems. Since high resolution imaging or spectroscopic fol¬ 
lowup is limited and expensive, we optimized the lens 
system redshift distribution to give maximal cosmology 
leverage. The optimization code under fixed resources 
such as spectroscopic time (e.g. to measure the lens 
galaxy velocity dispersion to constrain the lens mass 
model) is computationally fast and efficient, with its algo¬ 
rithm generally applicable to many astrophysical studies 
and figures of merit. 

The sculpted distribution delivers a nearly 40% im¬ 
provement in dark energy figure of merit, and a factor 
two tighter constraint on the Hubble constant, than a 
uniform redshift distribution. Low redshift systems are 
found to be particularly preferred, and there is no need 
to spend followup time on lenses with z > 0.5. 

Systematics enter as correlated quantities within a 
given lensing system, and as model systematics common 
to many systems. We examined both in an illustrative 
model that captures key aspects. A systematic uncer¬ 
tainty floor somewhat spreads out the optimal redshift 
distribution, but preserves the advantage of low redshift. 
We then demonstrated the effects of both a systematic 
explicitly evolving in redshift and one caused by popu¬ 
lation drift between different lens system types, such as 
double image vs quad image systems. If the systematic 
level for one population is larger than for the other, we 
can quantify by how much the followup resources are bet¬ 
ter spent on the more accurate population. 

As wide field surveys deliver 1000-10000 strong lensing 
systems, the issue of followup will become a key limita¬ 
tion, and these optimization tools can significantly im¬ 
prove the cosmological leverage. Similarly, as our mea¬ 
surement of strong lensing systems improves, the illustra¬ 
tive systematic correlation model here will become more 
realistic and enable more sophisticated trade studies re¬ 
garding low vs high redshift, or double vs quad image 
systems. This will further optimize future surveys to use 
strong lensing time delay distances as a unique cosmo¬ 
logical probe. 
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